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Abstract 

The stability of the magnetically charged Reissner-Nordstrom black hole solu- 
tion is investigated in the context of a theory with massive charged vector mesons. 
By exploiting the spherical symmetry of the problem, the linear perturbations 
about the Reissner-Nordstrom solution can be decomposed into modes of definite 
angular momentum J. For each value of J, unstable modes appear if the horizon 
radius is less than a critical value that depends on the vector meson gyromagnetic 
ratio g and the monopole magnetic charge q/e. It is shown that such a critical 
radius exists (except in the anomalous case q = with < 5^ < 2), provided only 
that the vector meson mass is not too close to the Planck mass. The value of the 
critical radius is determined numerically for a number of values of J. The instabil- 
ities found here imply the existence of stable solutions with nonzero vector fields 
("hair") outside the horizon; unless q = 1 and g > 0, these will not be spherically 
symmetric. 
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1. Introduction 



The Reissner-Nordstrom solution of the coupled Maxwell-Einstein equations 
describes a spherically symmetric black hole endowed with electric or magnetic 
charge. It may remain a solution when the Maxwell theory is embedded in a larger 
theory. In particular, a trivial extension of the Reissner-Nordstrom solution gives 
a magnetically charged solution to a spontaneously broken SU{2) gauge theory. 
However, it has been shown [1] that if the black hole horizon is sufficiently small 
(roughly, less than the Compton wavelength of the massive vector meson) this 
solution develops a classical instability. This instabihty suggests that there should 
be some other, stable solution with the same mass and magnetic charge. This 
solution has been found [2]; it is a gravitationally perturbed magnetic monopole 
with a black hole inside it. 

The spontaneously broken gauge theory can be viewed as a special case of a 
much wider class of theories obtained by varying the vector meson magnetic mo- 
ment and by allowing a less restrictive set of interactions. Although such theories 
will not in general be renormahzable, they are perfectly sensible classically. The 
magnetically charged solutions of these theories have been studied [3] both with 
and without gravity. Without gravity, the solutions are singular unless certain re- 
lations between the couplings are satisfied. When gravity is added, there are both 
Reissner-Nordstrom-type solutions and, for certain ranges of parameters, magnet- 
ically charged solutions with nontrivial matter fields outside the horizon. In this 
paper we analyze the stability of the Reissner-Nordstrom solutions against small 
perturbations. We determine the critical value of the horizon radius below which 
the Reissner-Nordstrom solution is unstable, and for which the less trivial new 
solutions are guaranteed to exist. In a future paper we will show how the results 
of the present analysis can be used to construct such new solutions. 

The Reissner-Nordstrom metric may be written in the form 

ds^ = -B(r)df + B-\r)dr^ + r^dO^ + sin^ Odcj)'^, (1.1) 
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where 



r / \ r 

with Qm and (5s being the magnetic and electric charges, respectively. We consider 
here the purely magnetic case, with Qe — ^ and Qm — (l/^-i where the Dirac 
quantization condition restricts q to be either an integer or a half-integer. Without 
loss of generality, we may assume that q is positive. In order that the singularity 
at r = be hidden within a horizon, we assume that M is greater than or equal 
to the extremal mass Mcxt = V^{q/ e)M-ph and hence that the outer horizon 
r+ = rjj > ^/47r{q/e)Mpy. The magnetic charge gives rise to a radial magnetic 
field 

F^^ = ^sin^. (1.3) 

All other components of F^,^ vanish. 

We want to consider the effect on this solution of introducing a massive charged 
vector field VF^. This field might be a non-Abelian gauge field which has acquired a 
mass through the Higgs mechanism (this was the only case considered in Ref. [1]), 
but this need not be the case. For our purposes here, we will only need the terms 
in the action which are quadratic in W. These take the form 



-\\Df,Wu-D,W^\^-rr?\W^,\' 



(1.4) 



Here F^jj = d^Ap — di,A^ is the electromagnetic field strength while = d^—ieAfj^ 
is the gauge covariant derivative. The W mass, m, might depend on the value of 
some scalar field 0, as in the Higgs case, but this dependence can be ignored here; m 
is then the value corresponding to the vacuum value of 0. The last term shown is a 
magnetic moment term, with the gyromagnetic factor g arbitrary. The integration 
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is over the region of spacetime outside the horizon of the black hole; it is this region 
that is relevant for investigation of black hole stability. We must therefore impose 
boundary conditions at the horizon, which will be discussed later. 

The Reissner-Nordstrom solution is trivially extended to this theory by setting 
— and setting all scalar fields at their vacuum values everywhere. We now 
want to consider fluctuations about this solution. These can be divided into several 
decoupled sets, comprising the metric and electromagnetic perturbations, the scalar 
field perturbations, and the vector field perturbations. The first set has been 
analyzed in detail [4] and shown not to lead to any instability. Because the scalar 
fields are assumed to be at their vacuum values, it is easy to show that small 
perturbations in these fields also cannot give rise to any instability; the proof 
parallels the proof [5] of the no-hair theorem for a free massive scalar field. The 
analysis of the remaining perturbations, those in W^, is the subject of this paper. 

The remainder of this paper is organized as follows. In Sec. 2, we exploit 
the spherical symmetry of the unperturbed solution by expanding Wf^ as a sum of 
eigenmodes of angular momentum and then decomposing the quadratic action into 
a sum of terms, each of which involves only modes of a given angular momentum. 
In Sec. 3, we show how both Wt and Wr can be eliminated, thus reducing the 
problem to one involving at most two radial functions for each value of the angular 
momentum. In Sec. 4 we show that instabilities develop if the horizon radius rjj is 
less than some critical value rcr, and obtain some analytic bounds on this critical 
value. Section 5 describes our numerical results for this critical radius. Some 
concluding remarks are contained in Sec. 6. 
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2. Angular momentum decomposition of modes 



The spherical symmetry of the underlying solution allows the perturbations to 
be decomposed into modes of definite angular momentum. Because of the non- 
trivial transformation properties of the monopole vector potential under rotations 
(or, more fundamentally, because of the extra angular momentum arising from any 
chargc-monopole pair) one cannot use the ordinary spherical harmonics for this 
decomposition, but must instead employ monopole spherical harmonics. We begin 
this section by reviewing the properties of these harmonics. Although our eventual 
interest is in a curved space-time, for the sake of compactness we adopt here a flat 
space notation with vectors denoted by bold face. 

For a scalar field of charge e, the appropriate spherical harmonics are the Wu- 
Yang monopole harmonics Yq^Mij^i'P) [6], which are eigenfunctions of the orbital 
angular momentum 

L = -ir X D - gf (2.1) 

obeying 

l?YqLM^L{L + l)YqLM 

(2.2) 

LzYqLM = MYqLM- 

The quantum numbers L and M take on the values 

L = q, q + 1, ... 

(2.3) 

M = -L, -L + 1, ...,L. 

These harmonics form a complete orthonormal set, with 

TT 27r 

J d^Y*^MYqL,M' = J smOde J d(t)Y*i^MYqL'M' ^ ^MM'- (2-4) 



For spin-one fields one must use monopole vector harmonics, which are eigen- 
functions of J = L -I- S. In general, there will be several of these corresponding to 
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the same values of and Jz- We will find it convenient to use a set [7], denoted 
by G^qjM, in which these are distinguished by the value of f • S. These obey 

JzC^Sm = M<^% (2-5) 
and are normalized so that 



^Or^(A)* ^(V) _ 5 J J' 5 MM' Sxx' , . 

^^''^qJM ■ ^qJ'M' " " 



The allowed values for the total angular momentum quantum number J occur 
in integer steps beginning with (? — 1, unless g = or 1/2, in which case the 
minimum value is J = g. As usual, the values of M run from — J to J in integer 
steps. Generically, there are three sets of harmonics for each J, with A = 1, 0, and 
— 1. However, some of these are absent if J is equal to or — 1. For J — q — 1 
there is only a single family of harmonics, with A = 1; these harmonics, which will 
be of particular importance in our analysis, have the important property that both 
their covariant curl and their covariant divergence vanish. For J = g > there are 
two families of harmonics, with A = 1 and 0. Finally, for J = g = (i.e., ordinary 
vector spherical harmonics), only A = is present. 

The last line of Eq. (2.5) shows that the A = harmonics must be purely radial; 
the normalization condition (2.6) then fixes them (up to an arbitrary phase) to be 

C^M = \i\jM. (2.7) 

Since, as one can easily show, vector harmonics with different values of A are 
orthogonal at each point in space, those with A = ±1 can only have angular 
components. 
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Wc will need the formulas for the covariant curls of the vector harmonics. For 
writing these, it is useful to define 



J2 = J( J + 1) - g2 



(2.8) 



and 

The curls can then be written as 

.(±1) _ 

.(0) 



D X C 



qJM 



±- 



(2.9) 



(2.10) 



qJM '^-^qJM 

(Note that = if J = g — 1, implying that the corresponding harmonics are 
indeed curl-free, as asserted above.) 

We will also make use of the identity 



qJM 



(2.11) 



Let us now proceed to the decomposition of the quadratic action. We begin 
with a mode expansion of VF^. Recalling that the A = harmonics are purely 
radial while those with A = ±1 are purely angular, and using Eq. (2.7), we write 



cx; J 

j=qM=-J 
^ oo J 



(2.12) 



J=q M=-J 
oo J 



J=q~l M=-J 



.(1) 
'qJM 



oo J 



J=q+1 M=-J 



where an index a represents either 9 or (p. We now insert these expansions into the 
action and keep only terms quadratic in W^. Using the properties of the harmonics 
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given above (and keeping track of the additional metric factors arising in our curved 
space-time) we obtain a sum of terms, each corresponding to definite values of J 
and M: 

oo J 

'S'quad = X/ X/ "^qu^d- (2-13) 



J=q-1 M=-J 



The action for the lowest angular momentum, J = g — 1, is particularly simple. 
As noted above, there is only a single multiplet of vector harmonics, and these 
have the properties of being purely angular and of having vanishing covariant curl. 
One finds that 



oo 



2 2| ^ |2 I 19 I r \2 



(2.14) 



rH 



(Here, and henceforth, we omit the superscripts JM on the various coefficient 
functions.) 

The general case is much more complicated, with 



qJM 

quad 



dt I dr{ |6-ra'P + 



+ 



B 



/+ - /c+a 



- B 



—m 



f+ - -k+h 



f- - -k-b 



/_ - k-a 



1 , ,2 

~ ;:2 1^+/+ ~ ^-/-l 



\U\' + \U\' + B\b\'-r'^\a\ 



+^o/.i^-i/-n}. 



(2.15) 

This simplifies a bit when J — q, since the terms containing /_ are then absent. 
As one would expect, it reduces to Eq. (2.14) if J is set equal to g — 1 and a, 6, 
and /_ are set equal to zero. 

Our task is now to analyze the behavior of perturbations governed by the 
various S^^. As a first step, we will now show that the functions a and b (corre- 
sponding to Wt and Wr, respectively) can be eliminated from the analysis. 
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3. Elimination of Wt and Wr 

The time derivatives of Wt do not enter the action (2.15). Hence, it (or a{r,t) 
in the reduced action (2.15)) is a nondynamical constrained field that can be elim- 
inated. We proceed as follows. If the dynamical fields /_, and b are assembled 
into a vector 

^+ ^-,bY, (3.1) 



the action (2.15) may be written in the form 



SiiL - I df l dr 



oo 

j dr zh + {z^Fa + a*F^i) + a*Ga - z^Hz , (3.2) 



where 



(3.4) 



and i7 is a 3 X 3 matrix, the exphcit form of which we do not need at the moment. 
(In going from Eq. (2.15) to Eq. (3.2) we have performed an integration by parts 
and dropped the surface terms. In the next section we will determine the boundary 
conditions on the fields; with these boundary conditions, the neglect of the surface 
terms is justified.) 

Variations with respect to a*{r,t) and z^r,t) yield the field equations 

= F^i + Ga (3.5) 

and 

0^z + Fd + Hz. (3.6) 
Using the first of these to solve for a, and then substituting the result into the 
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second equation, we obtain 

= (^7 - FG-^F^^ z + Hz. (3.7) 

Instabilities of the Reissner-Nordstrom solution correspond to solutions of 
Eq. (3.7) that have a time dependence of the form e"^ with real a > 0. If the 
operator acting on z is positive definite, these solutions will be in one-to-one cor- 
respondence with the negative eigenvalues oi H. To show that this operator is 
positive, we define the projection operator 

P = F{F^F)-^F^ (3.8) 

and write the operator in question in the manifestly positive form 

/ - FG-^F^ = {I - P) + PF 

(Note that, although the unstable modes correspond to the negative eigenvalues 
of if, the actual value of a is not given by the corresponding eigenvalue unless 
the unstable mode lies in the subspace spanned by 7 — P. This will not affect 
our calculations, because we are addressing only the question of the existence of 
instabilities.) 

Our problem is thus reduced to the study of the positivity of the potential 
energy 



oo 




(3.10) 



which is independent of Wt- Because the radial derivative of h does not enter the 
action, a further reduction is possible. We write as the sum of two integrals. 



1 



1 



Ftp. (3.9) 
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one of which contains only /_|_ and /_, and the other of which is a perfect square: 



I \2 , \tl\2 



i/;r+i/: 



1 



\k+f'+ + k-f: 



I |2 



oo 

+ J drB 



(3.11) 

The second integral on the right is clearly nonnegative. Furthermore, given any 
f+{r) and f~{r), it is always possible to choose b so that this integral vanishes. 
Hence, a necessary and sufficient condition for instability is that there be configu- 
rations of f+{r) and /-(r) for which the first integral is negative. We denote this 
quantity by E^^^, and write 



E. 



JM 
quad 



oo 

Jdr B f^Kjf + /t (^^I - ^) / 



(3.12) 



where / = (/+, /_)-^ and we have defined matrices 



Kj^I- 



1 



_|_ j2 



( k\ k+k-\ 



k+k- k' 



J > q 



(3.13) 



and 



Vj 



f-kl + f k^k- \ 



k+k- 



k 



2 19 



J > q. 



(3.14) 



For J — q — 1 OT q, matters are simplified by the absence of /_, so that Kj and 
Vj are numbers rather than matrices. Using the fact that kl. is equal to and q, 
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respectively, for these two cases, we obtain 



Kq-i = 1, = f (3-15) 



and 



^, = ^^. V;=(^. (3.16) 

(The results for J = g — 1 could, of course, be obtained directly from Eq. (2.14).) 

The positivity of Kj will be of importance later. This property is obvious for 
J — q — 1 and J — q, while for J > g it follows from the positivity of the two 
eigenvalues 

ki — 1 and k2 — — k ^tk. (3-17) 



4. Existence of instabilities 

Our task is now to analyze the positivity of the potential energies E^^^. In 
doing this, it convenient to transform to a tortoise coordinate x defined by 



dx 
dr 



B-\r). 



(4.1) 



Spatial infinity, r = oo, corresponds to a; = oo, while the horizon radius rjj 
corresponds to x — — oo. In terms of this coordinate 



E. 



JM 
quad 



— I dx 



dx dx 



(4.2) 



where 



Ueff = B{r) 



2r Vj 

m I - -J 



(4.3) 



Instabilities correspond to the existence of functions f{x) which for which the 
energy functional is negative. This is equivalent to the existence of a negative 
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energy bound state in the spectrum of the Hamiltonian 



(4.4) 



Because K is a positive matrix, a necessary condition for the occurrence of an 
instabihty is that U^ff be negative (or have a negative eigenvalue) in some region 
outside the horizon. This requires that Vj have a positive eigenvalue and that 
the horizon radius rj{ be less than a value ro(J) = a/m. For J = q — 1, this gives 
the requirements that g be positive and that rji be less than 



f 1^ 9Q -1 

'^o(?- 1) = W ■ 



(4.5) 



For J — q, instability can arise only if g > 2, while 



/ N i{9-'^)(l -1 



(4.6) 



U J > q the two eigenvalues of Vj are 



v± 



-J^±^J^ + g(g-2)q^ 



(4.7) 



If < (7 < 2, both of these eigenvalues are negative and no instabilities can arise. 
For values of g outside this range, is positive and 



ro{J) 



-J^ + ^J^ + g{g-2)q^ 



1/2 



m 



J > q. 



(4.8) 



Note that for fixed q and g, ro(J) is a decreasing function of J. 

Because of the gradient term in the energy, these conditions by themselves are 
not sufficient for the existence of an unstable mode. In fact, the horizon radius 
must be less than a critical value rcr(J) that is somewhat smaller than ro(J). In 
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the next section we will describe the numerical determination of rcr- Before doing 
so, we show that an unstable mode will always arise if rjj is sufficiently small, 
i.e., that there is in fact a positive rcr- (Our proof implicitly assumes that is 
greater than the extremal horizon radius. This is a significant constraint on the 
existence of rcr only if m is comparable to M-p\.) To see that this should not be 
immediately obvious, consider the flat space analogue of our problem in which B{r) 
and Kj{r) are both identically equal to unity, and the integration is restricted to 
a region r > r* > with the boundary condition that /(r*) = /(oo) = 0. Prom 
the identity 

^ 2 ,^|2" 



/*|/f^/ 



dr 



2r 



4^2 



(4.9) 



it is easy to see that the energy functional is positive, no matter how small r* is, 
if V (or its smallest eigenvalue) is less than 1/4. 

The situation is different in curved space, essentially because the metric factor 
suppresses the effects of the gradient terms near the horizon. To demonstrate this, 
let us consider the functional 



E' = dr 



B\ff+[m'-l]\f? 



(4.10) 



with y > being equal to the positive eigenvalue of Vj and 



(4.ii; 



For a given value of rjj, the Reissner-Nordstrom B{r) is less than B{r) everywhere 
outside the horizon. From this, together with the fact that Kj{r) < 1, we see that 
if the potential energy functional (4.2) is positive, then so must be E'. Conversely, 
if the Hamiltonian H' that we obtain from E' has a bound state, then so must the 
H of Eq. (4.4), and hence there will be an instability. Our aim is to show that 
such a bound state develops if rjj is made sufficiently small (but still positive). 
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with m and V held fixed. Dimensional analysis shows that this limit is equivalent 
to holding V and fixed and letting m become arbitrarily small. We will find it 
more convenient to use this latter method. 

Defining a tortoise coordinate y hy dy/dr — B~^{r), we have 



H' = -^ + U{r{y)), (4.12) 



with 



U = B{m^-^]. (4.13) 



Now write 



where 



U^Ui + U2 + Us, (4.14) 



Ui{y) = (4.15) 
10, r > Nrn 



{U{y), NrH<r<ro 
otherwise 



U2{y) = ^^'' , . (4.16) 



r 0, r < ro 
Us{y) = , ^ (4.17) 
{U{y), r > ro, 

with N being some fixed, large number and ro = VV /m. It is assumed that m is 
small enough that ro ^ Nth ^ ru- 

Consider the Hamiltonian H'q obtained by replacing U{y) by Ui{y). Because 
Ui{y) is everywhere negative on the real line, a standard result of elementary 
quantum mechanics guarantees the existence of at least one negative energy bound 
state. Let us denote the energy of this state by Eb = —c?, and let V'd/) be the 
corresponding normalized wave function. We do not need a detailed expression for 
but only the fact that it has a smooth nonzero limit as m — > 0. This can be 
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seen by noting that setting m = changes the value of Ui{y) at any point by a 
fractional amount of at most w?{Nr}^Y /V = {Nri^/rof' <^ 1. For a sufficiently 
large value of rg, the range r > tq will be entirely in the classically forbidden 
region. Hence, in this region, -0 will be of the form 'il^{y) — Ae~'^y with A only 
weakly dependent on m. 

An upper bound on the minimum eigenvalue £'min of H' can now be obtained 

by using ip as a trial wavefunction: 



oo 

Emin<EB+ J dy{U2 + U3M'' 



-oo 

oo 

|2 



y{ro) (4.18) 



oo 



<-Q;2 + m2 I dylAlh-^'^y 



2 

2a 



For sufficiently small m (i.e., large tq and y(ro)), the second term is smaller in 
magnitude than the first. Hence E'min is negative and we have proven our result. 



5. Numerical Results 

In the original coordinates, we have an energy functional (3.12) for the per- 
turbations. As we have shown, the black hole is unstable if and only if there exist 
functions /+ and /_ such that the appropriate functional is negative. We find 
such functions by looking for solutions to the eigenvalue equation with negative 
eigenvalue for the Hamiltonian defined by the functional.''^ The critical radius rcr is 
the value of the horizon radius above which there are no negative eigenvalues and 

1 This equation can also be found by varying the energy functional with respect to / and 
inserting a Lagrange multiplier to implement the normalization constraint. 
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below which there is at least one. When inspecting the equations, it is simpler to 
consider the case where r//, and thus the boundary condition, is fixed and m varies. 
(The relevant dimensionless parameter is mrH-) It is then clear that decreasing m 
strictly and smoothly decreases the energy, and thus the energy eigenvalues; the 
critical value of mru will be where the lowest eigenvalue crosses from positive to 
negative. 

A critical point in the parameter space, then, will be one for which the Hamilto- 
nian has (generically) one zero eigenvalue and no negative eigenvalues. In general, 
the wavefunction has two components, and /_. For J < q, however, there is 
only one component and so to find the critical point we search for values of the 
parameters such that the solution to the eigenvalue equation with eigenvalue zero 
has the following properties: 

1) The solution exponentially decays to zero at r — > oo. We can solve the 
equation analytically for large r; we choose the solution that is decaying. 

2) The solution goes to a finite value at r th- This is actually a fairly 
restrictive condition; a typical solution blows up at the horizon. That this is the 
correct boundary condition can be seen by looking in tortoise coordinates: In these 
coordinates a solution with a negative energy eigenvalue decays exponentially as 
X — > — oo. In the limit that the negative eigenvalue goes to zero, this exponential 
decay is arbitrarily slow. 

3) The solution does not change sign. 

Properties (1) and (2) are the appropriate boundary conditions; property (3) 
ensures that there are no negative eigenvalues. (In one dimension, the ground state 
wavefunction is the unique energy eigenfunction with no zeroes.) 

Prom the argument in the previous section, we see that there will always be a 
bound state for sufficiently small rjji or equivalently, sufficiently small m. Thus 
we can fix the other parameters and search for the critical value of m for which 
there is a solution to the differential equation that satisfies (1) through (3). We 
enforce condition (1) by integrating inward from r ^ ro, choosing initial values 
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to approximate the decaying exponential we expect at large r. (The precise ini- 
tial values chosen are not important, since the decaying exponential solution will 
quickly dominate as we integrate to smaller r.) If the solution changes sign before 
we reach rjj, we know that m is below the critical value; if the solution diverges 

2 

as r ^ rjj without changing sign, we know that m is above the critical value. 
We can use this information to quite rapidly and accurately zero in on the critical 
value of m (or, equivalently, ru). 

For J = q — 1, the energy functional is given by 



E. 



(q- 
quad 



00 

1)M_ f 



dr 



m 



(5.1) 



thus the (zero eigenvalue) equation to be solved is 

d 



(5.2) 



In this case there is only one other dimensionless parameter, gq. (We make the 
approximation r_ = 0, good for m <^ M-p\.) We have calculated the critical mass 
for a variety of values of this parameter. 

For J — q, we have 



quad 



dr 



B 



2 

1-1)1 

r'^m? + q 



TH 



which gives the zero eigenvalue equation 



d_ 
dr 



B 



r'^vr? 
r'^rr? + q 



- 0) 



+ m 



2r2 



{g - 2)q 
2r2 



1/4 



(5.3) 



(5.4) 



In this equation q appears separately from g; in principle, the critical mass depends 
on l)oth. 



2 The algorithm we used in fact searches for a sign change, but if it does not find one, it 
tests the value of df /d\n(r — rn) very close to rn- If this is positive, it is assumed that the 
solution will in fact change sign closer to the horizon than the integrator can proceed, and 
thus m is still below the critical value. 



17 



In the case J > q, the wavefunction has two components. This gives a two- 
component equation, which does not have a simple relation between nodes and 
ground states and which would require tuning of the relative weights of the two 
components, as well as the mass. In this case, we are able to calculate bounds 
on the energy by minimizing one-component functionals; these energy bounds give 
bounds on the critical parameter. 

The energy functional for J > q is given by Eq. (3.12). Consider a basis in 
which Vj is diagonal, and where the upper left element is v+, the only positive 
eigenvalue of Vj. (See Eq. (4.7)) The matrix Kj has two positive eigenvalues, ki 
and k2 (Eq. (3.17)); ki is the larger of the two eigenvalues. 

Now, we are searching for the minimum energy configuration for the func- 
tional E[f] (subject to some normalization constraint). We call this configuration 
/o. Now let /i be a configuration which minimizes the functional, subject to the 
additional constraint that the lower component of / is zero (in this particular 
basis). Clearly E[fi\ > -Bf/o]- Now consider the modified functional E-[f] ob- 
tained by substituting the smaller eigenvalue of Kj, k2 (actually k2l), for Kj. 
Clearly E-[f] < E[f] for all /, so if we minimize E- with a configuration /2, 
E-[f2] < E[fo]. Thus we have bounded E[fo] above and below. But both E[fi] 
and E-[f2] can be calculated with a single component functional, which allows 
us to put bounds on the critical parameters without tuning the two-component 
equation. 

Now if we write 

then in the basis where Vj is diagonal (with the positive eigenvalue v+ at the upper 



3 Since the upper left component of Vj is the only one which is positive, the energy will never 
be negative if we require that the upper component of / is zero. 
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left), 
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(5.6) 



where v — ^/ J^^ + q^g{g — 2). So, to calculate the minimum energy with the lower 
component constrained to be zero, E[fi\, we minimize the functional 
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rH 



where C„„ — ^ + ^[^^ + ^^id ~ 2)1 upper left component of the matrix 

C defined above. If we tune m to solve the zero eigenvalue equation associated 
with this functional, we obtain a lower bound on the critical mass. We bound the 
energy from below (and thus the critical mass from above) with the described 
above. Since both Vj and Kj are diagonal in this functional, the upper and lower 
components of / completely decouple: E_ = E_[fu] + E_[fi], where /„, // are the 
upper and lower components of / respectively. Because E-[fi] is positive definite, 
we can set fi — and just use the functional 
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JM 
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rH 



Solving the zero eigenvalue equation associated with this functional gives an upper 
bound on the critical mass. 

We have calculated and graphed some critical values (or bounds thereon) of 
mrui where m is the W mass and rji is the outer horizon radius, for varying g, 
q, and J. (We have calculated these assuming that r_ — 0, i.e. a Schwarzschild 
solution, which will be approximately correct as long as m is far below the Planck 
scale. Provided that the black hole is not very near extremahty, the calculated 



19 



values will not be significantly affected by a non-zero r_. Handling the extremal 
case is considerably more delicate, and we do not discuss it here.) 

First note that the upper and lower bounds on rcr calculated for J > are 
quite close, which justifies calculating only these bounds rather than doing the 
full, two-component calculation. All the curves show the qualitative dependence of 
the critical radius on g (and J). Most importantly, it appears from our numerical 
results that r^r is a decreasing function of J; recall that ro( J) has the same property. 
Thus, as we reduce the radius of a black hole with fixed charge, it will be the mode 
with the lowest angular momentum that first becomes unstable. If g is positive 
and g > ^, then this will always be the mode with J = q — 1. If g is negative 
the J < q modes never become unstable (see Eqs. (4.5) and (4.6)) and it is the 
J — q + 1 mode that first becomes unstable. (We do not display results for < 0. 
For the J > q modes, making the substitution g ^ 2 — g does not change the 
calculated upper bound on rcr, although it does change the lower bound. The 
qualitative behavior is quite similar.) Finally, if g = |, then there is no instability 
if < < 2; for 5( > 2, it is the J — q — ^ mode that is the first to become 
unstable. 

It is useful to compare the behavior if rcr with that of ro. In Figure 2 we plot 
the ratio of these quantities. Note that for J > q, we find that rcr ~ tq for all 
values of g, with the ratio varying from about ^ to 1. For J — q — 1 we find 
similar behavior for large g, but for very small g, we find that the relationship 
between g and rcr becomes linear, i.e. that rcr oc Vq. (Indeed, for small g, we have 
iTir-cj- = {mro)'^, with unit proportionality.) We have no analytical explanation for 
this behavior. Note also that the values of Tcy/tq for J = 1 arc very close to the 
values for J = 0, when g > 2. Although the differences between these curves are 
very small, they are much greater than the estimated errors in our calculation. 
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6. Conclusion 



We have investigated the hnear stabihty of magnetically charged black holes 
in a theory with a massive charged vector field. The Reissner-Nordstrom solution 
is stable against small perturbations of the electromagnetic, gravitational, and 
scalar fields, so linear perturbations of only the W field need be considered. We 
expanded the W field in angular momentum eigenstates, the monopole vector 
harmonics. For each value of J, we found the critical horizon radius at which 
an unstable mode appears. This radius (in units of the W mass) depends upon 
the angular momentum of the mode, J, the charge of the black hole, q, and the 
magnetic moment of the W field, g. This makes quantitative the arguments of 
Ref. [3] concerning the stabihty of Reissner-Nordstrom black holes in this theory. 

There are some unanswered questions in this analysis. In particular, when 
the critical horizon radius is close to the radius of an extremal Reissner-Nordstrom 
black hole with the same charge, these calculations must be modified. Indeed, if the 
critical radius we have calculated here is less than the extremal radius, there may 
be no instability at all. Because of the quite different nature of the coordinates near 
the extremal horizon, our algorithm for calculating the critical value was ineffective, 
and so this case was not directly addressed. Another outstanding question is the 
origin of the linear behavior of the critical radius for J = q — 1 near g = 0. 

In general, it is the mode with the lowest J that becomes unstable at the largest 
th- Consider the evaporation of a Reissner-Nordstrom black hole with magnetic 
charge q/e. Once the horizon radius falls below rcr for the lowest mode, we expect 
that the black hole will be described by a nontrivial classically stable solution with 
given and q, rather than the unstable Reissner-Nordstrom solution. If g = 1, 
then the lowest mode has J = and is spherically symmetric, so we expect that 
this new stable solution will also be spherically symmetric. Suppose instead that 
q ^ 1. The lowest mode will have J > and the new classically stable solution will 
not be spherically symmetric. The further evolution of the black hole depenmds 
on whether or not the theory possesses finite energy magnetic monopoles that are 
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not black holes. If it docs, then as the horizon contracts with further evaporation 
we expect that the black hole will decay by emitting monopoles and settling down 
to a g = 1 spherically symmetric solution. Eventually the horizon will disappear, 
leaving behind a simple monopole. If there are no such monopole solutions, the 
endpoint of the decay will be an asymmetric extremal black hole. 
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FIGURE CAPTIONS 

1) The critical value of mrH as a function of for ^ = 1 and J — 0, 1, 2, and 5. 
For J — 2 and 5, upper and lower bounds are plotted for mr cr, rather than 
the actual values. 

2) The ratio of Vcr to rg as a function of g for the values of J shown in Figure 1. 
As in Figure 1, bounds are plotted for J > g = 1. Note that the curve for 
J = 1 is largely obscured by the curve for J = (for g > 2). 
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